Thermal Stability of Strained Nanowires 
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Stranski-Krastanov strained islands undergo a shape anisotropy transition as they grow in size, 
finally evolving toward nanowires. This effect has been explained until now via simple energetic 
models that neglect thermodynamics. We investigate theoretically the stability of strained nanowires 
under thermal fluctuations of the long side. We find phase transitions from nanowires back to 
nanoislands as the temperature is increased and as the height of the nanostructure is raised or 
lowered and we predict regions of phase coexistence. Our results are general, but explain recent 
data on the growth of erbium silicide on a vicinal Si surface. 

PACS numbers: 61.46.-w, 65.80.+n, 68.35.Rh, 68.55.ag 



The investigation of two-dimensional nanoislands has 
provided fundamental understanding of the mechanism 
of epitaxial growth. Instabilities in heterostructure 
growth are exploited to form self-organized nanodots and 
nanowires that are of much smaller size than similar 
nanostructures realized on lithographically pre-patterned 
substrates [TJ Typically, competition of plastic and 
elastic relaxation processes in strained epitaxial layers 
drives their formation [3j 0j. In Stranski-Krastanov 
growth [3J [5], where interactions between adatoms and 
of adatoms with the surface are comparable, both dis- 
located or single crystalline islands can grow after the 
initial wetting layer has reached a critical thickness. For- 
mation of single crystalline rather than dislocated islands 
depends on the material, the amplitude of the lattice mis- 
fit and the possible interplay of the Asaro-Tiller-Grinficld 
instability [51 [7]. It can be controlled in many ways, e.g. 
by undulations and/or patterning in the surface layers of 
the substrate [8] or by tuning areal island densities and 
temperature [9]. 

It has been predicted and observed experimentally 
that dislocation free islands in early stages of growth 
can undergo a shape anisotropy transition as their size 
becomes critical [TU] [HI H2]- Eventually their growth 
becomes quasi-one dimensional and they form perfect 
nano-objects with well defined height and width, that 
is, nanowires. Previous theoretical treatments have been 
based on simple energetic models of the island/ wire that 
include competition between strain relief and surface en- 
ergy. Although useful, these models completely neglect 
thermal fluctuations which are important in such quasi- 
one dimensional systems because high entropic gains of 
the (thermodynamically extensive) long sides can lead to 
various instabilities. 

Our objective is thus to introduce thermal effects in 
the theory within the context of a statistical mechanical 
model. We find that above a critical temperature the 
wire does indeed become unstable and this temperature 
depends on the threshold of the plastic relaxation of the 
film. Moreover, we find that the wires become unstable 



as the deposition continues at constant temperature. We 
explain the observations of a recent set of experiments 
on the growth of erbium silicide nanostructures on the 
vicinal Si(001) surface that was studied at different post- 
annealing temperatures and coverages [131 HH HH] • 

We first review the Tersoff and Tromp model of the 
shape transition [TO] on which more detailed studies have 
been based [UJ . By considering both surface energy 
and strain relaxation Ref [10] arrives at the following ex- 
pression for the energy for unit area of a rectangular is- 
land of dimensions L, X 
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the parameters d, e Q depend upon surface energy and 
strain. The superficial density of energy in (JT|) is mini- 
mized by a square island of size X = L = d and energy 
dre a . Yet, depending on the growth condition, islands 
can grow metastably larger than the optimal size d 2 . 
Minimizing at fixed area A — XL, Tersoff and Tromp ob- 
tained a critical area, A c = (ec?) 2 , below which symmetric 
shapes, or X = L = \/~A, are energetically favorable and 
above which a symmetry breaking occurs. That is, as 
the surface A increases beyond the critical value, one of 
the lateral sizes, say X, shrinks down from ed asymptot- 
ically approaching the equilibrium value d. The other 
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FIG. 1: The nano-wire with lateral edges xi(l), X2Q) allowed 
to fluctuate in the coordinate I, < I < L; x(l) = x\{l) — X2(l) 
is the relevant degree of freedom for our problem. Amplitude 
of fluctuations is largely exaggerated in figure. 
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side keeps growing eventually linearly as L ~ A/d. In- 
deed this asymptotic behavior can be extracted from ([I]) 
by taking the limit L — > oo and then obtaining the su- 
perficial energy for lateral growth 



V(x) = 



d xe 
-£o-ln— : 
x d 



(2) 



which has as expected minimum in x = d (we will from 
now on use lower case x to denote the finite edge, and 
upper case L for the long edge - "long" as in the thermo- 
dynamic limit, Fig. [TJ . 

We will explore the thermal stability of a long wire 
whose lateral size d minimizes the asymptotic energy in 
([2]) under fluctuations of its long edges. We allow the 
lateral edges x\(l), X2U) to fluctuate in the coordinate 
< I < L. The problem is simplified by the introduc- 
tion of x = xi — x 2 , x + = xi + x 2 - The elastic energy 
of those fluctuations is thus simply T = kx' 2 + kx' 2 2 = 
\kx' 2 + \kx' 2 . As the term in x + only contributes to 
equipartition, it can be disregarded. The partition func- 
tion in the continuum limit is a gaussian path integral 
over the allowed fluctuations 



Z((3) = []l[dx{l)]e-Ptf*[$°rm+*rMm] 
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Standard transformations |16[ I17j allow us to write it as 
the trace 
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of the one-particle, temperature dependent hamiltonian 
given by 



(5) 



The coefficient in front of the trace simply adds to the 
free energy the equipartitive contribution of free oscil- 
lators (A is the linear density of oscillators A = K/k, 
with K the microscopic spring constant). We thus only 
need to concern ourselves with the diagonalization of the 
(temperature dependent) effective Hamiltonian. In the 
limit of a long wire, or L/d — > 00, the density oper- 
ator in Q projects on the lowest bound state - pro- 
vided that Hp has one. One can reproduce the pre- 
vious calculation for the partition function and demon- 
strate that the probability of a lateral size x is given by 
p(x) — linij^ + \ip^ ' J (x)\ 2 where t/j { p J is the lowest bound 

state of the operator Hp + Jx, where J is a chemical po- 
tential. Clearly when Hp has no bound states - which 
might happen for (3 below a certain j3 c - then p(x) is 
flat, and we conclude that the elongated structure is not 
stable at that temperature. Thus the search for stability 
of long wires is reduced to the problem of existence of a 
bound state in a suitable 1-particle hamiltonian opera- 
tor, a situation reminiscent of localization/delocalization 
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FIG. 2: From top to bottom: the potential of |2| (solid line) 
with plastic limit set at x p = 8d (dashed line). Behavior of 
T c (in units of T = tt^ 1 \/Wke^d?) as a function of the plastic 
limit x n 



of one-dimensional interfaces of the solid on solid model 

nans iag. 

Before exploring this transition, one can immediately 
show a relationship between the critical exponent of x = 
J x\ipp(x)\ 2 dx ~ t~~P and that of the specific heat C ~ 
r a when t = (T — T c )/T c -> 0". Indeed, from Q we 
have F ~ TLE°/d, while for T - T~ one has E° ~ 



' . since it is the exponential tail of the bound 



t 2 ? 



or 



state around infinity that dictates x. Thus / 
a = 2(1 - 0). We show at the end that (3 = 1. 

The potential in ^ converges to zero slower than 1 jx 
at infinity while at zero it goes to infinity slower than the 
centrifugal energy; thus Hp, for any value of coupling 
constant (i.e. any temperature) has at least as many 
bound states (excluding S waves) as the hydrogen atom: 
infinite. Such a cursory look might convince the reader 
that nanowire formation should always be stable with 
respect to thermal fluctuations. In fact it is the long 
range behavior of the potential that commands possible 
instability, but ^ is clearly unphysical for large x; sim- 
ple intuition suggests that for large islands, energy den- 
sity must be constant, as plastic relaxation of the clas- 
tic strain becomes unavoidable when the islands grow 
in size. Relaxation via dislocations has been observed 
experimentally in nanoislands and has already found to 
promote many complex phenomena .21], only recently 
explored via simplified approaches |221 123] . Fortunately 
our scenario is simpler: islands grow in a single crystalline 
manner until they reach a threshold size in which dislo- 
cations appear. Once dislocations start to form, they can 
coalesce, move and make other dislocations more likely. 
It is thus both reasonable and simple to introduce an ad 
hoc plastic limit x p beyond which the energy of the island 
simply grows linear with size. 
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FIG. 3: Phase diagram (T, h) from Q and the experimental 
findings of Zhou et al. [TJj. Boxes denote nanowires, empty 
boxes regions of coexistence while diamonds correspond to 
nanoislands (the relative coverage is 0.29, 0.73, 0.87, 1.16, 
1.45, 2.03, 2.32, 2.90 Mono Layer or ML). We have chosen 
our parameters arbitrarily (yet reasonably) to provide a de- 
cent superposition with Cai's results: — 6 nm, i = 1 nm, 
T c ,mm = 548 °C. Note that for low h, where T c diverges, the 
region of coexistence of islands and wires is still experimen- 
tally accessible, as discussed theoretically in the text. Note 
also the region of nanowire stability with no transition for 
T < T c , min and for large heights h > 3.5 nm. 

We substitute the potential in ([2|, V — > V p , with 
V p (x) = V(x) when x < x p and V p {x) — V{x p ) = V p 
for x > x p . V p (x) is now a well for x > (Fig. [2]) 
and thus might or might not admit a bound state, de- 
pending on the value of (3. The Schrodingcr problem for 
the potential in ^ has not been solved yet, neverthe- 
less the critical temperature corresponding to the disap- 
pearance of the lower bound eigenstate can be obtained 
via a WKB approximation, which works well when T c is 
not too large [23] : 

T c = l —^dPg(^f). (6) 
For £ > 1, g(£) is defined by 

<?(£) = J ^lnM-iln(e^) ds (7) 

fo is the lower zero of the integrand. 

Equations ([6]) , Q link critical temperature to the plas- 
tic limit. Figure [2] depicts the behavior of T as a func- 
tion of the threshold for dislocation x p assuming all the 
other parameters are constant. One would deduce that 
in general stability can be achieved at higher temper- 
atures by pushing higher the single crystalline thresh- 
old. Yet the phase transition was obtained in the limit 
d/L I 0. Consider now Q, (|5| for L/d large but finite. 
For T > T c , where T c is the critical temperature, there 
are clearly no nano- wires, as there are still no bound 
states. When T < T c , the operator in ^ has (at least) 




FIG. 4: A series of STM images (all 300 X 300 nm 2 ) taken on 
the vicinal Si(001) surface with different Er coverages. The 
samples are post-annealed at 650 C for 10 min. (a) 0.29 ML. 
(b) 0.73 ML.(c) 2.03 ML. (d) 2.90 ML. From ref [T3]. 

one bound state, and yet the density operator does not 
just project on the ground state: it leaves some proba- 
bility p oc e~^ E "t>\ of a jump from the bound state to 
the continuum spectrum, a jump activated by the small 
"pseudo-temperature" d/L. In the proximity of T c , we 
obtain (see later) E% = ir 2 e d 2 (T c - T) 2 /2T C 3 ; p is thus a 
gaussian in T — T c of relative width 

AT C 1 \~TT 
l c 7r V e Ld 

which implies a gray zone of coexistence of both islands 
and wires of the order of 10 _1 T C (if we assume k of the 
same order as e d, L/d ~ 10 2 and x p /d ~ 10 — 100). 
Notice that the width in ^ increases with T c . That 
makes it possible for a region of coexistence of islands and 
wires to be accessible even when T c is beyond practical 
reach. This seems to be the case in ref [T3], as explained 
later. 

It is interesting to explore the behavior of the criti- 
cal temperature with deposition. In general the plas- 
tic limit as well as e , k and d depend on island height 
h, surface energy and strain. Tersoff and Tromp found 
d/e = 4>hcxp(T /2ch) and e Q — 2ch/d, where T is propor- 
tional to the surface energy, c to the strain energy and 
4> is a number which depends on the contact angle |10j . 
Clearly, the linear clastic constant k — xh is also pro- 
portional to the island height. In the simpler approach, 
the energy balance x p (y{x p ) — V(x p /2)) = T p provides 
a criterion for the formation of a defect of energy (per 
unit length) T p , which gives x p /d = 4exp (T p /2ch — 1). 
So, if we neglect the dependance of T p upon temperature 
we can relate critical temperature with the island height 
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and thus obtain the critical line 
2e</> 



T c (h) 



Xc h e 2ch g -e 2 <= h 



(9) 



In Fig. [3] we plot Q as a phase diagram for a 
nanowire/nanoisland under the above assumptions, along 
with experimental results of Zhou et al. [14 . One sees 
that below a certain temperature T Cim ; n the growth of a 
nanowire by further deposition is always stable. Above 
Tc,min, deposition will cause growth in height until a crit- 
ical value is encountered, above which the wire becomes 
unstable. Equation ^ also predicts a region of wire sta- 
bility for higher values of h. As Fig. [3]shows, the interval 
of instability in h decreases in size with the temperature, 
until it disappears below T C)m i n . Depending on materi- 
als, kinetic parameters and other physical circumstances, 
only a portion of the phase diagram of Fig. [3] might be 
accessible. 

Our predictions explain the experimental results of 
Cai's group on the self organization of erbium silicide on 
Si (001). They obtained phase transitions from nanowires 
to nanoislands, both when increasing post-annealing 
temperatures at fixed coverage [T5] and when increas- 
ing coverage (and thus the height of nanostructures) at 
fixed post-annealing temperatures as in Figs. [3j|4] [13 14J. 
In a yet different approach, the same group was able to 
control the efficiency of the chemical synthesis of the sili- 
cide via different deposition techniques, and could obtain 
nanostructures of different heights (2.0, 1.5 , 0.5 nm) at 
the same temperature: again transition from nanowires 
to nanoislands was observed as the height increased |15j . 
as well as the coexistence zone. 

Finally, let us tie some loose ends. The WKB method 
seems appropriate to compute the critical temperature 
but not for critical exponents. It is clearly so for a case 
that can be solved exactly, the square well. And yet, as 
the length scale diverges at the transition, one expects 
the actual shape of the well to become irrelevant to the 
exponents. With the approximation g(£) ~ (£ — l)/2, 
which works fine for £ of the order of 3-10, one finds 
that our critical temperature corresponds to that of a 
square well of depth e Q and size a such that 4e Q a 2 = 
£o(x p — d) 2 . The problem for the square well potential 
can be solved easily and returns for the ground state 
E° p = ^-e a d 2 (Tc ~ 3 T)2 , and then X/x~t and thus / ~ i 2 , 
for t < while / — for t > 0: / has no kinks, hence 
entropy is continuous at the transition and there is no 
latent heat, but a discontinuity in the specific heat. One 
can check that on the contrary a WKB study of a square 
well (and also for our truncated well) returns different 
exponents, X/x ~ t 1 / 2 , and thus a latent heat. 

In conclusion, we have studied the stability of nanowire 
fabrication under thermal fluctuations. We find phase 
transitions from wires to islands, both by increasing tem- 
perature and the height of the nanostructure. Plasticity 
of extended films plays a crucial role in our analysis. Our 



results show that the transition can be more or less sharp 
depending on the average length of the wires, and is pre- 
ceded by regions of coexistence of wires and islands. We 
thank Dr. Cai and her group for sharing figures and data. 
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